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Abstract 

We consider the energy and the polar angle distributions of J/ip''s produced 
via the color-singlet and the color-octet (^S'^^'*) mechanisms in — > J/ip-^-X 
at LEP. Since both distributions of the J/V' produced via color-octet mecha- 
nism are significantly different from those via color-singlet mechanism, these 
observables can be used as tests of color-octet production mechanism for heavy 
quarkonia. We also discuss Z° — > T + X and W J/Tp{oiT) -|- X in brief. 
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Heavy quarkonium productions in high energy processes have long been studied in per- 
turbative QCD along with the color-singlet model. However, the recent measurements of 
J / if) and if)' productions at the Tevatron put serious questions on such theoretical approaches 
|1|]. In short, the data for — > ^z^' + X is larger than the color-singlet model by a factor 
of ~ 30. In order to resolve this discrepancy, Braaten and Fleming suggested a color-octet 
fragmentation of an energetic gluon into if)' by emission of soft gluons [@]. This approach 
requires one parameter, (OlOg (^S'i)|0), the matrix element of a color-octet, dimension- six, 
four-quark operator in the Nonrelativistic QCD, which is the effective field theory for QCD 
relevant to the heavy quarkonium physics 0. Braaten and Fleming showed that if one 
fixes this parameter to fit the measured cross section, then the observed spectrum in is 
reproduced. Then, Cho and Leibovich [|] extended this approach in order to go beyond 
the fragmentation picture (which is applicable only to the limit p\ ^ "^^O^ including other 
color-octet intermediate states such as ^Sq~^^ and ^-Pj=o,i,2*^^'* in terms of two more NRQCD 
matrix elements, (O|0g (^5*0)10) and (OlOg ('^Po)|0). This is especially relevant for the case 
of inclusive T{nS) productions at the Tevatron, since the pt in this case is not that large 
compared to mx and thus the fragmentation picture is not valid any more. Overall de- 
scriptions for ip^'^ and T(nS) productions at the Tevatron seem reasonably good by now, 
although some color-octet matrix elements are shown to lead too large contributions to J/ ip 
photopro duct ions [|] and B ^ J/ip + X 0]. 

In view of this, it is important to test the idea of color-octet mechanism in the heavy 
quarkonium productions in places other than at the Tevatron. One can think of the 5— wave 
charmonium production in B decays [^,0, the angular dependence of J/ip production cross 
section in the e~^e~ annihilations ||], J/ip + X PJTO[|, J/ip productions at the 7p 

collision [§,0, the associated J/?/) + 7 hadroproduction ||rT| and many others. All of these 
have formed an active research field recently . 



In this letter, we reconsider the inclusive heavy quarkonium production in decays, 
which is accessible at the current LEP experiments. This study is partly motivated by 
the recent report by OPAL collaboration that they have observed excess of events for 



^0 _^ T{nS) + X for n = 1,2, 3, larger than theoretical expectation by a factor of ~ 10, 
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J2 B{Z' ^ T(n5) + X)exp ^ (1.2^°:^ ± 0.2) x 10-^ (1) 



n=l 



compared to the 6— quark fragmentation contribution 



B{Z^ T{nS) + X)frag ^ 1.6 x 10"^ (2) 



n=l 



Similar excess were also observed in the prompt J/ip and ip' production in Z° decays, 
although the experimental errors are quite large |jl5| : 



^ J/^ + X) = (3.0±2.3) X 10-^ (2.6x10-^) (3) 
^7/;' +X) = (2.2±1.5) X 10"^ (1.1 x 10"^) (4) 



where the results for the heavy quark fragmentations [11^ are shown in the parentheses 
for comparison. The issues of the branching ratios and the energy distributions of J/ip in 

^ J /ip + X have been addressed already in Refs. ^ and . In this letter, we suggest 
another observable, the polar angle distribution of J /ip (relative to the e+e~ beam directions 
at LEP) for the color-singlet and the color-octet contributions, in order to provide another 
independent check of the idea of color-octet mechanism in heavy quarkonium productions. 

In order to study the angular distribution of J/ip produced in the Z^ decays at LEP, it 
is convenient to define S{E) and a{E) as 

S,s(£,)(l + COS^S.), (5) 

following the case of e+e" ^ 7* ^ J/V- + X at CLEO energy, ^ 10.6 GeV [|,|Tg. The 
subscripts 1,8 denote the color-singlet and the color-octet contributions, respectively. Since 
we are concerned with the polar angle (cos^^^) distribution of J/il) in the rest frame of Z^ 
as well as the energy distribution, we have to include the polarization of the Z^ produced 
through the e^e~ annihilations at LEP (with ^/s = Mz)- This effect can be conveniently 
described by replacement of e^{Z)e*^{Z) in the decay process by the density matrix p^^ 1 



where is the 4-momentum of . 

The vector polarization and the tensor polarization Q}^^ in Eq. (|^) can be obtained 
in the case of t^{k2) er{k\) — > Z^ process as following [|1^] : 



V 



2- 



1 / ,., Z^Z^X A^A" 



where A'^ = {k\ — k'l)^. 

After obtaining the invariant amplitudes for the color-singlet and the color-octet pro- 
cesses, we have performed the symbolic manipulations over the squared amplitude with 
the above polarization tensors and integrated over the phase space, using REDUCE and 
Mathematica. The resulting expressions for Sx^'^iE) and a\^'^{E) are rather lengthy, and will 
be shown separately elsewhere [T^. Our formalism with the density matrices for polarized 



particles can be applied to other processes where Z^'s are produced via different mechanism 
from the e^e~ — * For such processes, the density matrix for Z° will be different from 
Eqs. (7). The explicit relation between the density matrices of J/^/^ and Z^ in Z^ J/tp + X 
for an arbitrary Z^ polarization will be discussed there, too. 

Let us first consider the conventional color-sing let contribution to Z° ^ J/ip+X at LEP. 
The most dominant contribution at LEP energy scale is known 0,[1O| to come from Z^ — >■ 
{cc){'^S[^^) + C + C, which can be further approximated into the (anti)c— quark fragmentation 



into J/ip as |T^. In this work, we employ the full amplitude for the color-singlet J/^p 
production without using the fragmentation approximation. The relevant Feynman diagrams 
are shown in Fig. 1. The corresponding amplitude can be obtained as usual in the covariant 
gauge. For simplicity, we do not show the amplitude for the color-singlet J/ip production 
in Z^ —y J/ip + X explicitly. However, the corresponding expressions for 5*1 {E) and ai (E) 



will be given later ||18| 



Next, we consider Z^ qq + g followed by (7 — > '^Si^'^ J/ip+{soii hadrons) with 
q = u,d, c, s, b. The relevant Feynman diagrams are shown in Fig. 2. Other diagrams are 
higher orders in f^, or suppressed by the short distance factor mg/Mf, and thus irrelevant 



to Z'^ ^ qq + J /ip( OTip' , T(nS)). It is convenient to use the following variables : 

s^iZ-Pf = {q, + q2f, (8) 
t^{Z-q,f = iP + q2Y, (9) 
u = {Z-q2f = {P + qiy, (10) 

with s + t + u = M| + + 2m^, and P, qi and q2 are the 4-momenta of J/ip, q and q, 
respectively. The color-octet contribution of J/ ip production via the process Zq qq + J/ip 
factorizes into short and long distance pieces as 

M{Zo ggJ/^)octct 

= M{Zo qqg*{^ CCp^f ^])),hort distance X M{ccf S[^^] J/ V^),o„g distance. (H) 

In terms of s,t,u variables, the amplitude for the process Zq — > qqg*{—* ccfs[^^]) is given 
by 

M{Zo ^ qqg*{-^ ccp^f ^Doctet = ^| ^(-Mf + -MfXl^X^P), (12) 

Mj/^ sm dw cos 6w 

Mr = u{<li)TaY^^^^^^r{9v + 9Al,)v{q2\ 
u — 



MT = u{q,)T^r{9v + 9Al,)\^^Yv{q,). (13) 



"q 



where gy and gA for the up-type and the down-type quarks are given as 

(+i ~ ' forg = M, c, t, 

{gv^gA) = < ' ' (14) 

[ (~i + I ^^^^ +i) ' for g = d, s, b. 

The long distance contribution due to the color-octet fragmentation of a gluon into the 

color-singlet heavy quarkonium J/ip, ip' or T(nS) is described in terms of a parameter 

A4s{nS), which depends on the heavy quarkonium state one considers. Cho and Leibovich 

have determined these parameters for nS = T{nS), J/ijj and ijj' @JT^ : for example, 
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MsiTinS)) = J2 \M{bb[^S[''^] ^ r{nS))\^ = 6.4 x 10"=^ GeV^ (15) 



n=l n=l 



MsiJ/^ilS)) = \M{cc[^sl^^] J/^)\^ = 0.68 X lO'^ GeVl (16) 



After one averages and sums the squared amplitude over the initial and the final spin states 
using Eqs. (6), and performs the phase space integration, one easily finds various distribu- 
tions as well as the branching ratios. 

We first show the z{= 2E^/Mz) dependence of ai^s{E^) for the color-singlet (in dashed 
curve) and the color-octet (in solid curve) case respectively in Fig. 3 (a). We note that 
the ai{E) abruptly changes its sign from ~ —0.8 to -|-1 near the end point of 2; ~ 1. This 
originates from the smallness of 5 = 2M^/Mz ~ 6.6 x 10^^. ^ In Fig. 3 (b). we show the 
a{Er) for the case of T(IS') + X. The behavior near z ^ 1 becomes less abrupt here, 

as a result of the not-too-small 6 ~ 0.42 in this case. We have checked that this abrupt 
behavior near z ^ 1 dies out as we increase 6. 

Our energy distributions shown in Fig. 4 (a) agree with those obtained by |^ and [jlO 



which is a check of our calculations based on the polarized decays. The color-octet 
mechanism (shown in the solid curve) produces softer J/ip's compared to the color-singlet 
mechanism (shown in the dashed curve). Most J/ip^s are predicted to have energy around 
~ 5.5 GeV, if the color-octet mechanism works. Similar is true of the T energy distribution 
shown in Fig. 4 (b). 

The cos 6^ distributions(6'^ is the angle between the initial electron and the final J/i/j 
direction) can be obtained by integrating (5) over the heavy quarkonium energy E^. The 
color-singlet contribution predicts the polar angle distribution to be 



d cos O^p 
dVi 



d cos Of 

On the other hand, the color-octet mechanism predicts 

° — /'OK c 1 n-2 



(5.4 X 10-^) (1 + 0.92 cos^ B) MeV. (17) 
(9.7 X 10-^) (1 + 0.73 cos^ Q) MeV. (18) 



d cos 



(25.8 X 10"^) (1 + 0.34 cos^ Q) MeV. (19) 



^ For the case of the J/V' production at CLEO energy through e+e — > 7* J /il) + X [Q] |16|, 
the corresponding 5 = M^/E = 3.1 GeV/5.29 GeV ^ 0.59 is not small. 



7-p 

(36.7 X 10"^) (1 + 0.27 cos^^) MeV. (20) 



d cos 9r 

Adding these two contributions together, we get 



(31.2 X IQ-"^) (1 + 0.44 cos^ 9) MeV. (21) 



d cos 9^ 

- (46.4 X 10-3) (1 + 0.37 cos^^) MeV. (22) 



d cos 9r 

Therefore, we find that the color-octet mechanism makes the angular distributions fiatter 
compared to that by the color-singlet mechanism alone. The polar angle distributions for 
Z ^ J/ip + X and Z — > T + X are shown in Figs. 5 (a) and (b). For each plot, the color- 
singlet, the color-octet and the sum of the two are shown in the dotted, the dashed and the 
solid curves, respectively. The difference in the 9^ distributions between the singlet and the 
octet contributions are not so prominent compared with the energy (E^) distributions shown 
in Figs. 4 (a) and (b). However, using the number of ~ 600 J/ip^s from Z'^ ^ J/ip + X a.t 
LEP, it would be possible to distinguish (21) from (17) from the cos6'^ distribution |T^ (at 
least for the J/ip production in Z° decays). Thus, the measurement of cos 9^ distribution 
could constitute another independent test of the idea of color-octet mechanism as a possible 
solution to ijj' anomaly at the Tevatron. Since the measurement of the polarization of ijj' is 
not easy, our option may be more viable experimentally. 

Finally, we can obtain the branching ratios for Z^ J/iIj{ot T(IS')) + X through 
the color-singlet and the color-octet mechanisms by integrating (17) -(23) over the whole 
9tp. The results are shown in Table 1. The color-octet contribution is comparable to that 
through the color-singlet contributions ( or equivalently, the c(or 6)— quark fragmentation 
into J/iIj{ot T(IS')) considered by Braaten, Cheung and Yuan [|14|)- When summed over the 
light fiavor quantum numbers (g = u, d, c, s, b), the former actually becomes larger than the 
latter by a factor of 4 — 10. For the ip' , T(2S) and T(35') production rates, the color-octet 
contributions should be scaled down by a factor of 0.3, 0.6 and 0.14, respectively, coming 
from A4s{2S)/M.8{lS). Our results on the branching ratio agree with the results obtained 
in Refs. [§] and [jlO|, which is another check of our calculations based on polarized Z° decays. 



Our analysis can be easily extended to the W decays, W ^ qq + g followed by (7 — >■ 
J/ip, (or ijj' ,T(nS)). Using the unitarity of the CKM matrix element V^^', one can easily 
show that 

B{W ^ J/^/j + X)octet ~ 2 X 10-^ (23) 



compared to the charm quark fragmentation in the singlet model [14 



B{W J/^ + X)frag ~ 4 X 10"^ (24) 

Again, the color-octet mechanism dominates the heavy quark fragmentation in the color- 
singlet model in W decays, and this may be observed at LEP200. For Upsilon production, 
we have 

B{W T{1S) + X)o,tet ~ 3.7 X 10^1 (25) 

In conclusion, we have considered the color-octet mechanism in the heavy quarkonium 
production in decays. Both T+X and J/tp+X are enhanced due to the color-octet 

mechanism by a factor of ~ 10 or so. Also, the energy and the polar angle distributions of 
J/ip's produced via color-octet mechanism are drastically different from those (via charm- 
quark fragmentation) in the color-singlet model. These two distributions can provide key 
tests for the idea of the color-octet gluon fragmentation which was invented to solve the ^' 
anomaly at the Tevatron. This idea might be the simultaneous solution to the ip' anomaly 
at the Tevatron and the excess of events in the channel Z° ^ T + X at LEP, if the observed 
energy and polar angle distributions follow predictions made here and by other works 
T0| . Similar phenomena occur in the case of inclusive W decays into J/ip^ip' ,T(nS), and 



the excess of W T + X would constitute another test of the color-octet mechanism. 

A remark is in order before closing. Some recent works on the color-octet contributions to 
J/ip photoproduction and B J/ip + X show that some of the color-octet matrix 
elements (via the ^S'o''^^ and ^Pj^^^ ) determined from the J/ip productions at the Tevatron |^ 
might have been overestimated by an order of magnitude. However, in the J/ ip productions 
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in Z° decays considered in this work, the color-octet ^Si ' channel is most dominant, with 
other color-octet channels being essentially negligible. Therefore, our predictions would be 
relatively stable. 
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FIGURE CAPTIONS 



Fig.l Feynman diagrams for the color-singlet mechanism for — > (cc)(^S'} ) + cc. 

Fig. 2 Feynman diagrams for the color-octet mechanism for qq -\- J /il) with 

q — u,d, c, s, b. 

Fig.3 The ai,s{z) of the (a) J/^) and (b) T as functions of z = 2Eh/Mz {H = ip, T) 
for the color-singlet (in dashed curve) and the color-octet (in solid curve) case respectively. 

Fig.4 The energy distributions of the (a) J/V' and (b) T. The color-octet contribution 
and the color-singlet contributions are shown in the solid and the dashed curves, respec- 
tively. 

Fig. 5 The polar angle distributions of the (a) J/-^ and (b) T (relative to the e+e" beam 
direction at the Z— peak). The color-octet contribution, the color-singlet contributions and 

the sum of the two are shown in the dashed, the dotted and the solid curves, respectively. 
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TABLES 

TABLE I. Branching ratios for ^ qq + J/ijj{oT T(IS')) for q = u,d,c,s,b with = 5 
MeV, 171(1 = 10 MeV, iris = 150 MeV, rric = 1.5 GeV and = 4.9 GeV ,and for the long-range 
matrix elements J2l=i Mi{T (nS)) = 2.3 x IQ-^GeV^ J2l=i MsiT (nS)) = 6.4 x IQ-^GeV^, 
Mi{J/ip{lS)) = 6.1 X 10-2GeV2 and MsiJ^lS)) = 6.8 x IQ-^GeV^. 



Q 


tf[Z qq + J/ip) 


±>[A qq+ L [lb)) 


u 


4.1 X 10~^ 


5.6 X 10"^ 


d, s 


5.4 X 10-^ 


7.2 X 10-6 


c 


3.9 X 10-5 


5.6 X 10-6 


b 


4.4 X 10-5 


6.5 X 10-6 


octet sum 


23.2 X 10-5 


32.1 X 10-6 


singlet 


5.6 X 10-5 


9.6 X 10-6 


(fragmentation) 


(6.5 X 10-5) 


(14.2 X 10-6) 


total 


28.8 X 10-5 


41.7 X 10-6 
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FIGURES 




FIG. 1. Feynman diagrams for the color-singlet mechanism for (cc)(^5'| ) + cc. 




FIG. 2. Feynman diagrams for the color-octet mechanism for — > qq + J/ip with 
q = u,d, c, s, b. 
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FIG. 4. The energy distributions of the (a) J/V' and (b) T. The color-octet contribution and 
the color-singlet contributions are shown in the solid and the dashed curves, respectively. 
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FIG. 5. The polar angle distributions of the (a) J/ if) and (b) T (relative to the e^e" beam 
direction at the .Z— peak). The color-octet contribution, the color-singlet contributions and the 
sum of the two are shown in the dashed, the dotted and the solid curves, respectively. 
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